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Accurate  knowledge  of  heat  transfer  is  required  as  it  determines  the 
structural  requirements  of  high  speed  vehicles.  Thermal  design  of  the  struc- 
tures and  material  selection  depend  heavily  on  engineering  estimates  of  the 
local  heat  transfer  rates.  In  the  case  of  turbulent  boundary  layers,  one  has 
to  rely  heavily  on  experimental  data,  since  no  theories  are  able  to  predict 
the  turbulent  heat  transfer  rate  without  incorporating  a large  amount  of 
experimental  data.  A recent  evaluation  of  such  empirical  theories  for  pre- 
dicting heat  transfer  in  constant  pressure  turbulent  boundary  layers  at  super- 
sonic and  hypersonic  Mach  numbers  is  given  by  Hopkins  and  Inouye  in  Ref.  1. 

These  empirical  theories  are  not  completely  adequate,  however,  to  pre- 
dict the  heating  rates  on  actual  geometries,  since  extensive  viscous-inviscid 
interactions  will  be  present  and  cause  heating  rates  many  factors  higher  then 
the  associated  flat  plate  values.  These  high  local  heating  rates  can  pose  a 
major  design  difficulty  as  discussed  by  Korkegi  in  Ref.  2. 

Many  of  these  viscous-inviscid  interactions  are  caused  by  shock  waves 
interacting  with  the  turbulent  boundary  layer  on  the  body  and  result  in  signi- 
ficantly higher  heating  rates.  Many  experimental  studies  have  been  conducted 
to  determine  the  magnitude  of  these  heating  peaks,  compared  to  the  flat  plate 
values,  plus  their  location  with  respect  to  the  provoking  agency.  Although 
many  of  these  studies  have  been  carried  out  with  nominal  two-dimensional  flow 
configurations,  there  have  also  been  a number  of  studies,  e.g.  Neumann  and 
Token,  Ref.  3,  that  do  consider  a particular  kind  of  shock  wave/boundary  layer 
interaction  which  has  essentially  three  spatial  dimensions.  Several  empirical 
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correlations  of  heat  transfer  amplification  factors,  with  respect  to  flat 
plate  values,  have  been  reported  for  various  viscous-inviscid  interactions. 
These,  however,  are  generally  limited  to  particular  geometries  and  test  condi- 
tions because  they  have  not  incorporated  a sufficiently  deep  physical  under- 
standing of  the  phenomena  involved.  Such  hit  and  miss  empiricisms  almost 
require  every  new  geometry  or  flow  configuration  to  be  experimentally  tested 
before  the  interaction  effects  can  be  appraised  properly.  This,  however,  can 
be  a very  costly  and  time  consuming  proposition. 

Recently  there  has  been  a series  of  developments  in  the  general  area 
of  numerical,  computational  techniques  and  high  speed  computers  pointing  the 
way  to  solving  the  full  time-averaged  Navier-Stokes  equations  on  a routine 
basis.  Fundamental  problems  remain  since  these  computer  codes  incorporate 
empirical  information  in  the  form  of  the  turbulent  Reynolds  stresses,  the  so- 
called  turbulence  models.  Thus  it  would  appear  that  these  elaborate  computa- 
tional techniques  suffer  the  same  kind  of  deficiencies  as  the  more  simplistic 
empiricisms,  e.g.,  for  predicting  peak  heating  rates.  These  drawbacks  can 
hopefully  be  overcome  based  upon  the  argument  that  the  high  level  of  sophisti- 
cation will  give  the  turbulence  models  a universal  applicability  depending 
only  upon  a basic  set  of  non-dimensional  parameters.  It  is  in  this  area  of 
computational  fluid  dynamics  that  large  bodies  of  sufficiently  accurate  ex- 
perimental data  for  well  defined  problems  are  needed  to  guide  the  development 
of  these  turbulence  models. 

The  objective  of  the  present  study  is  two-fold.  The  first  one  is  to 
provide  a detailed  set  of  surface  data,  complemented  by  a complete  survey  of 
the  associated  flow  field,  from  which  detailed  flow  models  can  be  constructed 
which  may  lead  to  a better  understanding  of  the  physical  phenomena.  This 
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would  be  particularly  useful  for  obtaining  direct  design  estimates  of  the 
particular  problem  under  consideration.  A second  and  more  far-sighted  objec- 
tive is  to  describe  and  document  a well  defined,  physically  significant 
problem,  that  has  been  investigated  experimentally,  which  can  guide  the 
development  and  validate  future  computer  codes  and  turbulence  models. 

With  these  objectives  in  mind,  the  present  study  considers  the  inter- 
action flow  field  and  surface  phenomena  created  by  the  three-dimensional, 
experimental  configuration  depicted  in  Fig.  1.  The  coordinate  system  used 
for  the  present  problem  is  shown  as  Fig.  2.  This  study  is  actually  a con- 
tinuation of  the  earlier  work  by  Oskam,  Bogdonoff  and  V as.  Ref.  4,  the  main 
difference  being  that  the  region  under  study  is  extended  into  the  streanwise 
corner  formed  by  the  shock  generator  and  the  tunnel  test  wall.  The  instru- 
mentation rows  have  therefore  been  taken  parallel  to  the  Y direction  instead 
of  being  parallel  to  the  X direction. 

This  document  is  basically  divided  into  three  parts.  The  first  one 
describes  the  environment  under  which  the  tests  have  been  conducted.  Special 
attention  is  paid  to  the  definitions  of  the  physical  boundary  conditions  under 
which  the  experimental  data  are  obtained.  This  is  of  particular  importance 
for  the  heat  transfer  studies  where  results  have  been  obtained  for  a uniform 
temperature  test  wall  on  which  the  variable  heat  transfer  rates  have  been 
measured.  A second  part  deals  with  description  of  the  apparatus  used  in  ob- 
taining the  various  data.  The  process  of  data  reduction  and  the  estimates 
of  data  accuracy  are  presented  in  the  third  part.  Much  attention  has  been 
given  in  this  respect  to  possible  limitations  of  the  techniques  employed  and 
the  resulting  accuracies  are  tabulated. 
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The  final  data  has  been  plotted  in  10  detailed  appendices,  but  can 
also  be  obtained  from  detailed  tabulations  as  given  by  Oskam,  Vas  and 
Bogdonoff  in  Ref.  5. 
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2.  EXPERIMENTAL  PROGRAM 

2 . 1 General  Object! ves 

In  an  attempt  to  conduct  an  experimental  program  which  gave  the  optimum 
amount  of  relevant  data,  special  attention  has  been  paid  to  an  exhaustive 
series  of  details  which,  hopefully,  give  the  final  results  a more  easily  under- 
stood meaning  and  eventually  will  lead  to  a more  basic  understanding  of  the 
physical  phenomena  involved. 

A special  wind  tunnel  test  section,  designed  and  constructed  during  the 
earlier  work  by  Oskam,  Bogdonoff  and  V as  (Ref.  4),  was  shown  to  be  an  effi- 
cient experimental  configuration  which  has  a relatively  thick,  high  Reynolds 
number,  turbulent  boundary  layer,  of  superior  transverse  uniformity  compared 
to  previous  studies  performed  with  tunnel  wall  boundary  layers  as  discussed  in 
Section  2.5.  The  maximum  shock  generator  angles  were  limited,  however,  by  the 
interaction  resulting  from  reflection  of  the  oblique  shock  wave  at  the  tunnel 
sidewall  causing  the  tunnel  to  stall  at  shock  generator  angles  larger  than  14°. 
This  test  section  features  two  large  observation/instrumentation  ports  (12  inch 
diameter)  which  rendered  flexibility  for  mounting  optical  windows,  surface 
instrumentation  and  probe  drives  covering  the  complete  area  of  interest  with- 
out affecting  the  shock  generator  model  configuration. 

The  region  studied  in  the  course  of  this  investigation  extends  from  the 
area  under  the  shock  wave  towards  the  axial  corner  formed  by  the  test  surface 
and  the  shock  generator,  including  the  turbulent  boundary  layer  on  the  shock 
generator  itself.  The  practical  importance  of  this  axial  corner  interaction 
in  high  speed  flow  stems  from  the  observation  of  high  local  heating  rates  in 
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the  immediate  vicinity  of  the  shock  generator.  A new  probe  drive  mechanism 
has  therefore  been  constructed,  accompanied  with  angle  sensitive  probes, 
total  temperature  probes  and  a static  pressure  probe,  and  has  been  employed, 
to  resolve  the  flow  field,  including  all  three  velocity  components,  up  to 
the  surface  of  the  shock  generator. 

Special  attention  is  given  in  this  study  to  check  the  validity  of  the 
experimental  techniques  employed.  The  heat  transfer  technique  is  described 
in  detail  together  with  the  precise  definition  of  the  coefficients  employed 
to  present  the  data. 

The  question  of  data  accuracy  has  extensively  been  discussed  in  this 
document  since  the  final  analysis  of  the  data  will  heavily  depend  on  the  accu- 
racy that  has  been  achieved.  The  accuracies  of  the  various  types  of  data  are 
given  in  terms  of  estimated  uncertainties  which  should  be  considered  to  be 
upper  bounds  since  they  were  obtained  by  considering  the  largest  possible  errors 
which  could  only  occur  at  certain  locations  of  the  region  studied. 

2.2  Test  Environment 

2.2.1.  Windtunnel  facility 

The  Aeronautical  Engineering  Department  at  Princeton  University  started 
its  study  of  viscous  phenomena  in  supersonic  flows  in  1946.  This  interest  led 
to  the  design,  construction  and  operation  of  the  4x8  inch  variable  density 
supersonic  tunnel.  That  tunnel  had  specifically  been  designed  to  cover  the 
range  of  Reynolds  numbers  up  to  full  scale  and  provided  experimental  informa- 
tion in  a test  range  not  covered  by  any  other  supersonic  testing  facility  at 
that  time  (Ref.  6).  After  many  projects  had  successfully  been  completed  the 
facility  was  moth-balled  due  to  the  laboratory's  emphasis  on  higher  Mach  number 
studies. 
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More  recently,  however,  there  has  been  renewed  interest  in  high  Reynolds 
number  studies  because  of  difficulties  encountered  with  extrapolating  wind 
tunnel  data  to  actual  flight  conditions.  As  a result  the  original  high  Reynolds 
number  facility  was  put  back  into  operation  in  1970  but  with  a new  Mach  3 noz- 
zle and  longer  test  section  (Ref.  7).  In  this  new  108  inch  long  boundary  layer 
channel,  with  a cross  sectional  area  of  8 x 8 inches,  one  can  obtain  Reynolds 

numbers,  based  on  the  turbulent  boundary  layer  running  length,  of  the  order  of 

q 

10  at  a tunnel  stagnation  pressure  of  500  psia.  These  newly  constructed  test 
sections  are  flexible,  easily  modified  and  readily  adapted  to  testing  of  vari- 
ous configurations. 

The  high  pressure  air  for  the  facility  is  provided  by  four  Worthington 
four  stage  compressors,  driven  by  100  horsepower  induction  motors,  which  can 
pump  from  atmospheric  up  to  3200  psia.  A drying  system,  operating  at  3200  psia, 
removes  all  the  moisture  from  the  high  pressure  air  to  eliminate  the  possibility 
of  the  adverse  effects  of  condensation  shocks  in  the  supersonic  tunnel.  Four 
tanks,  situated  above  the  ground,  provide  a storage  capacity  of  2000  cubic  feet 
at  a maximum  pressure  level  of  3200  psia. 

The  stagnation  pressure  is  regulated  by  an  8 inch  hydraulically  actuated 
valve  during  tunnel  operation.  A control  system,  which  can  be  operated  either 
manually  or  automatically,  is  used  to  maintain  the  desired  stagnation  pressure. 

A tunnel  run  has  to  be  terminated  if  the  tank  pressure  drops  so  low  that  the 
main  regulat  ' valve  reaches  maximum  opening  so  that  the  desired  stagnation 
pressure  can  no  longer  be  maintained. 

Since  the  tunnel  is  not  provided  with  a heat  exchanger,  it  is  to  be 
expected  that  the  stagnation  temperature  will  change  from  day-to-day  as  well 
as  during  tunnel  operation.  The  day-to-day  change  is  determined  by  the 
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variation  in  the  ambient  temperature  since  this  controls  the  initial  tempera- 
ture of  the  storage  tanks.  The  actual  stagnation  temperature  level  in  the 
test  section  is  also  determined  by  the  Joule-Thompson  drop  across  the  main 
regulator  valve,  and  a significant  amount  of  heat  transfer  to  the  air  from 
the  inlet  piping,  in  addition  to  the  initial  tank  temperature.  This  leads  to 
a drop  in  tunnel  stagnation  temperature  with  time  during  operation.  The  obser- 
vations that  the  tunnel  walls  are  initially  at  room  temperature,  and  that  the 
stagnation  temperature  is  always  below  room  temperature,  explain  the  experienced 
heat  transfer  from  the  walls  to  the  gas  and  result  in  a drop  of  the  tunnel  wall 
temperature  with  time.  The  magnitude  of  the  stagnation  temperature  decrease 
being  slightly  larger  than  the  wall  temperature  decrease  with  time,  imposes  a 
permanent  driving  potential  for  the  wall  heat  transfer  during  operation.  This 
enables  one  to  study  heat  transfer  rates  on  tunnel-wall  mounted  models  without 
being  in  direct  control  of  either  the  stagnation  or  the  tunnel  wall  temperature. 

The  above  mentioned  variations  in  temperature  conditions  from  run  to  run 
and  also  during  a single  run  do  require,  however,  special  data  reduction  pro- 
cedures which  will  be  discussed  later.  It  should  also  be  noted  that  the  free 
stream  Reynolds  number  is  a direct  function  of  the  stagnation  temperature  and 
is  thus  subject  to  similar  variations. 

2.2.2.  Test  Section  and  Shock  Generator  Apparatus 

The  model  configuration,  which  was  used  previously  (Ref.  4),  consists  of 
a sharp  flat  plate,  8 inches  by  10  inches  long,  called  the  shock  generator, 
mounted  vertically  between  tunnel  top  and  bottom  wall  and  pivoted  about  a center 
close  to  its  leading  edge.  This  shock  generator  and  its  bearings,  which  are 
held  in  the  top  and  bottom  walls,  were  constructed  out  of  one  piece  to  avoid 
leakage  beneath  the  leading  edge  along  the  tunnel  wall.  The  remaining  edges 
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of  the  shock  generator  were  lined  with  nylon  to  provide  the  seal. 

The  main  support  of  this  generator  is  provided  by  an  arm  which  extends 
through  the  tunnel  side-wall  and  is  linked  to  a screw  mechanism.  The  housing 
for  this  mechanism  is  sealed  to  the  side-wall  of  the  tunnel  and  has  an  external 
drive  wheel,  so  one  can  manually  rotate  the  shock  generator  from  outside  the 
tunnel.  A range  of  shock  generator  deflection  angles  from  0°  up  to  14°  has 
been  used  during  the  investigation.  Since  sizeable  two-dimensional  models  at 
large  deflection  angle  are  particularly  susceptible  to  tunnel  starting  problems, 
the  starting  procedure  adopted  was  to  initially  set  the  shock  generator  at  an 
angle  not  exceeding  8°  and  to  increase  its  deflection  angle  after  the  supersonic 
flow  was  established  if  desired.  Permanent  tunnel  stall  was  experienced,  how- 
ever, above  14°. 

Four  static  pressure  taps,  located  in  the  surface  of  the  shock  generator, 
are  connected  with  flexible  tubing  to  a plug  in  the  mechanism  housing.  Moreover, 
two  thermocouples  are  spot  welded  to  the  generator  to  measure  its  temperature. 

The  coordinates  (with  Oq  = 0;  of  the  static  pressure  taps  and  the  thermocouples 
are  given  in  Table  I. 


TABLE  I:  SHOCK  GENERATOR  INSTRUMENTATION  LOCATION 


Static  Pressure 
Tap  No. 

Thermocouple 

Number 

X 

Inches 

Z 

Inches 

1 

1 

2.0 

2.0 

2 

5.0 

2.0 

3 

2 

7.0 

2.0 

4 

2.0 

6.0 
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Upstream  of  the  shock  generator  leading  edge,  one  thermocouple  is  located  in 
the  tunnel  top  wall  to  indicate  the  wall  temperature. 

A photograph  of  the  test  section,  installed  between  the  nozzle  and  a 
second  test  section,  as  well  as  a line  drawing  of  the  shock  generator  actuating 
mechanism  are  presented  in  Figs.  3 and  4,  respectively. 

2.2.3.  Test  Conditions 

The  test  section  was  mounted  in  a position  immediately  downstream  of 
the  nozzle,  where  the  flow  uniformity  is  optimal,  and  the  starting  problems 
are  minimal.  The  nominal  free  stream  Mach  number  at  a location  upstream  of 
the  interaction  is  2.94.  The  turbulent  boundary  layer  on  the  tunnel  wall 
entering  the  interaction  is  approximately  0.55  inches  thick,  depending  on  the 
exact  location  under  consideration.  Only  part  of  the  heat  transfer  studies 
were  performed  over  a range  of  stagnation  pressures  from  75  to  200  psia  (see 
Section  2. 4.4.4).  All  of  the  remaining  tests  were  carried  out  with  a constant 
stagnation  pressure  of  100  psia  held  to  within  0.5%  during  tunnel  operation. 

At  this  100  psia  level,  a maximum  running  time  of  7 minutes  can  be  sustained, 
if  the  storage  tanks  are  at  a pressure  of  3000  psia  initially. 

Before  we  state  the  various  temperature  test  conditions,  it  is  most 
helpful  to  appraise  some  of  the  temperature-time  variations  mentioned  earlier 
in  Section  2.2.1.  Figure  5 shows  these  variations  versus  time  for  a run  con- 
sisting of  two  separate  tests,  15  minutes  apart,  and  each  approximately  180 
seconds  long,  at  a stagnation  pressure  of  100  psia.  One  may  observe  from  the 
data  that  after  an  initial  temperature  start-up  process  of  about  30  seconds, 
the  wall  temperature  and  stagnation  temperature  decrease  almost  linearly  with 
time.  Approximate  values  are: 
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The  shock  generator,  however,  due  to  its  large  surface  area  and  smaller  heat 
capacity  compared  to  the  3 inch  thick  aluminum  tunnel  wall,  decreases  much 
faster  in  temperature  than  the  tunnel  wall.  After  a time  of  about  two  minutes 
it  asymptotes  to  a level  close  to  an  adiabatic  wall  temperature  which  is  pro- 
portional to  the  tunnel  stagnation  temperature.  The  actual  proportionality 
factor  depends  somewhat  on  the  shock  generator  angle  under  consideration,  and 
is  measured  to  be  0.956  and  0.944  for  shock  generator  angles  of  10°  and  4° 
respecti vely . 

Since  the  test  times  typically  varied  between  1 and  3 minutes,  it  is 
necessary  to  quote  values  of  stagnation  and  wall  conditions  time-averaged  over 
the  actual  length  of  a test.  These  averages,  and  the  resulting  Reynolds  number 
per  inch,  are  tabulated  in  Table  II  for  the  measurements  of  surface  static  pres- 
sure and  flow  field  quantities.  The  variations  of  these  averages  resulted  from 
the  day-to-day  change  in  ambient  conditions: 


TABLE  II:  TEMPERATURE  TEST  CONDITIONS  FOR  SURFACE 

STATIC  PRESSURE  AND  FLOW  FIELD  MEASUREMENTS 


445  + 18°R 


503  + 17°R 


Reinch_ 


1.76  x 10°  + 


The  test  times  for  heat  transfer  studies  were  significantly  shorter. 
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with  a maximum  of  30  seconds.  Actual  data  was  only  taken  for  about  7 seconds, 
resulting  in  an  insignificant  change  in  temperature  conditions  during  the  data 
taking  process,  thus  eliminating  the  need  for  averaging.  The  variations  from 
run-to-run,  however,  are  still  experienced  as  indicated  in  Table  III. 

TABLE  III:  TEMPERATURE  TEST  CONDITIONS  FOR  HEAT  TRANSFER  TESTS 


Instrumentation  and  Techniques 
2.3.1.  Eccentric,  Rotatable  System 

The  test  section  features  two  12  inch  diameter  observation/instrumentation 
ports  (see  Fig.  4)  which  render  a unique  flexibility  for  mounting  optical  windows, 
surface  instrumentation,  and  probe  drives  for  covering  the  complete  area  of 
interest  without  affecting  the  shock  generator  positioning  and  actuating  mech- 
anism. A line  drawing  of  the  eccentric,  rotatable  system  which  can  accomodate 
surface  instrumentation  as  well  as  probe  drives  is  shown  in  Fig.  6.  It  consists 
of  a large  12  inch  diameter  plug  carrying  a smaller  7 inch  diameter  plug  2 inches 
off  center.  This  smaller  plug  has  a 2 inch  diameter  hole  located  again  2 inches 
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off  center.  The  actual  instrumented  plugs  are  located  in  this  2 inch  hole. 

The  various  plugs  were  fabricated  to  close  tolerances  to  get  a flush  surface 
inside  the  tunnel  and  sealed  with  0-rings.  The  rotatable  system  has  been  pro- 
vided with  scales  on  the  outside  so  one  can  reach  a given  location  by  rota- 
tion to  the  desired  angles. 

2.3.2.  Static  Pressure  Plug  and  Scanivalve 

The  2 inch  diameter  surface  static  pressure  plug  contains  45  static 
pressure  holes  distributed  in  3 rows  of  15  each  (see  Fig.  7).  The  typical 
spacing  of  taps  in  a row  is  0.100  inches.  The  tap  diameter  is  0.032  inches 
for  the  2 outer  rows  and  0.024  inches  for  the  middle  row.  Although  there  is, 
in  practice,  no  upper  limit  to  the  instrumentation  density  that  can  be  obtained 
with  the  eccentric,  rotatable  system,  the  static  pressure  plug  has  been  de- 
signed in  such  a way  that  a simple  180°  turn  of  the  instrumented  plug  relocates 
the  taps  at  positions  midway  between  their  location  before  the  180°  rotation. 

This  effectively  produced  an  instrumentation  density  of  one  tap  per  0.050  inches, 
or  10  taps  per  initial  boundary  layer  thickness. 

The  45  surface  pressures  were  sampled  by  a 48-port  Model  48J4  Scanivalve 
which  was  connected  to  the  taps  by  plastic  tubing.  The  pressures  were  con- 
verted to  electrical  signals  by  a Pace  transducer.  Although  the  volume  between 
the  Scanivalve  and  the  Pace  transducer  was  minimal  and  the  Scanivalve  case  was 
vented  to  a reservoir  held  at  the  undisturbed  tunnel  static  pressure,  the 
scanner  still  required  a minimum  of  0.4  seconds  per  reading  to  avoid  problems 
with  pressure-time  lag  at  the  transducer  when  the  Scanivalve  was  stepped  through 
a pressure  gradient. 

An  absolute  upper  bound  to  the  errors  caused  by  pressure-time  lag  in 
the  scanner  has  been  obtained  without  effort.  The  180°  turn  of  the  plug  to 
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obtain  the  higher  instrumentation  density,  while  keeping  the  same  plastic 
tubing  connections  and  the  same  Scani valve  rotation  direction,  reverses  the 
apparent  sign  of  the  pressure  gradient  as  observed  by  the  Scani valve.  The 
true  static  pressure  distribution  is  thus  bound  by  the  two  sets  of  data  ob- 
tained for  different  plug  orientations. 

Two  of  the  four  static  pressure  taps  on  the  shock  generator  were  also 
sampled  by  the  Scanivalve;  the  remaining  two  were  measured  by  two  separate 
Pace  transducers  connected  to  spare  channels  of  the  electronic  data  system. 

2.3.3.  Heat  Transfer  Measurement 
2 . 3 . 3 . 1 . Introduction 

The  experimental  method  used  to  determine  the  local  convective  heat 
transfer  rates  was  quasi -transient  and  calorimetric  in  nature.  The  temperature 
decrease  with  time  of  a piece  of  copper  of  known  mass  and  exposed  surface  area 
gave  the  heat  transfer  rate  directly.  The  driving  potential  for  this  heat 
transfer  is  determined  by  the  tunnel  operating  characteristics  as  discussed 
in  Section  2.2.1.  To  prevent  the  small  pieces  of  copper,  called  slug  calori- 
meters, from  reaching  the  adiabatic  wall  temperature  shortly  after  tunnel  start, 
a tiny  jet  of  hot  air  was  impinged  on  the  rear  side  of  the  slug  calorimeter  to 
control  its  temperature  while  the  tunnel  was  running. 

The  heat  transfer  data  was  derived  under  the  condition  of  uniform  wall 
temperature  along  the  test  surface.  To  satisfy  this  condition,  the  slug 
calorimeter  temperature  was  raised  by  the  jet  to  about  40°R  above  the  sur- 
rounding wall  temperature  while  the  tunnel  was  running.  At  that  time,  the 
heating  jet  was  cut  off  and  the  high  speed  data  acquisition  system  started. 

After  about  3 to  7 seconds,  the  slug  calorimeter  had  reached  a temperature  well 
below  the  wall  temperature  of  the  tunnel,  and  the  data  system  was  stopped.  The 
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temperature-time  gradient,  from  which  the  heat  transfer  rate  was  calculated, 
was  then  evaluated  at  the  time  that  the  slug  calorimeter  temperature  equaled 
the  surrounding  tunnel  wall  temperature,  so  that  non-uniformities  in  test  sur- 
face temperature  were  avoided. 

The  cycle  of  heat  transfer  measurement  can  be  repeated  during  a tunnel 
run.  The  restriction  in  the  number  of  cycles  taken  was  determined  by  the  data 
processing  rate  of  the  acquisition  system. 

2.3. 3.2.  Hardware 

A photograph  of  the  2 inch  diameter  heat  transfer  plug,  which  could  be 
located  in  the  eccentric,  rotatable  system,  is  shown  in  Fig.  8.  This  instru- 
mentation  is  a later  development  of  two  similar  heat  transfer  plugs  used  in 
the  earlier  work  (Ref.  4).  The  detailed  geometry  of  the  slug  calorimeter 
arrangement  and  the  hot  air  supply  tube  is  shown  in  a line  drawing  presented 
in  Fig.  9.  The  eight  slug  calorimeters  are  each  0.088  inch  in  diameter  and 
0.040  inch  thick  and  are  separated  from  the  2 inch  brass  plug  by  insulators 
made  of  linen  based  phenolic.  Alumel-chromel  thermocouple  wires,  0.005  inch 
in  diameter,  were  used  to  measure  the  temperature  of  the  slugs.  Two  tiny  holes 
of  0.006  inch  diameter  were  drilled  through  each  slug  near  its  center  to  accom- 
modate the  wires.  They  were  soldered  to  form  a thermocouple  junction  within 
the  copper  slug. 

Small  tubes,  with  an  inside  diameter  of  0.040  inches,  are  located  within 
0.12  inches  of  the  back  of  the  slugs.  These  eiqht  tubes  were  scaled  up-and 
connected  to  a single  half  inch  diameter  supply  tube  within  the  2 inch  plug 
itself.  Compressed  air  was  tapped  from  the  shop  air  system  at  a maximum  pres- 
sure of  120  psia,  passed  through  a boiling  water  heat  exchanger  and  supplied 
to  the  half-inch  tube.  A three-way  valve  located  in  the  line  immediately  before 
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the  2 inch  plug,  could  exhaust  the  hot  air  to  the  tunnel  room  or  to  the  instru- 
mentation. 
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Two  wall  temperature  sensors  are  located  in  the  2 inch  brass  plug  be- 
tween the  8 slug  calorimeters.  All  thermocouple  wires  were  led  to  a single 
reference  temperature  bath.  Both  legs  of  each  thermocouple  junction,  the 
alumel  and  chromel  wires  were  connected  to  copper  lead  wire,  and  these  two 
new  junctions  (chromel -copper,  alumel -copper)  were  held  at  the  water-ice  point. 
These  copper  wires  were  connected  to  the  inputs  of  very  high  impedance  linear 
amplifiers.  This  procedure  of  referencing  both  legs  avoids  the  generation  of 
thermal  emfs  at  the  input  terminals  of  the  amplifier. 

2. 3. 3. 3.  Test  Procedures 

Periodic  checks  of  the  thermocouples  were  made  at  room  temperature. 

This  was  done  by  removing  the  instrumentation  from  the  tunnel  and  placing  it 
in  an  isothermal  enclosure  together  with  a mercury  thermometer.  At  this  known 
temperature,  the  output  of  the  data  system,  after  calibration  with  an  accurate 
millivolt  source,  was  compared  with  the  standard  for  alumel -chromel  thermo- 
couples as  determined  by  the  U.  S.  National  Bureau  of  Standards. 

The  actual  procedure  for  taking  heat  transfer  measurements  consisted 
of  starting  the  tunnel  and  keeping  the  slug  calorimeters  at  a temperature 
close  to  the  tunnel  wall  temperature.  After  the  shock  generator  was  rotated 
to  the  desired  angle,  and  the  slugs  were  heated  to  about  40°R  above  the  wall 
temperature,  the  data  acquisition  system  was  started  and  the  three  way  valve 
in  the  hot  air  supply  line  was  turned  to  exhaust  to  the  tunnel  room.  This 
suddenly  stopped  the  heating  process  and  consequently  the  slugs  decreased  in 
temperature.  After  about  7 seconds  the  data  system  was  stopped.  A typical 
slug  temperature-time  trace  during  this  process  is  shown  in  Fig.  10. 
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The  actual  speed  of  the  multiplexer  was  200  readings  per  second.  With 
10  channels  this  results  in  20  readings  per  thermocouple  per  second.  Although 


the  in-house  IBM  System  7 computer  can  record  as  many  as  1000  readings  per 
second,  and  store  up  to  1600  readings,  the  overall  data  through-put  rate  was 
limited  by  the  connection  between  the  System  7 and  the  main  computer,  which 
operates  on  a time  sharing  basis.  This  resulted  in  a tine-averaged  data  flow 
rate  to  the  main  computer  of  about  10  readings  per  second.  Thus,  short  heat 
transfer  runs  could  be  made  successively  with  time  intervals  of  about  5 minutes. 

The  spacing  between  the  slug  calorimeters  was  0.45  inches.  Their  loca- 
tion was  such  that  a simple  180°  turn  of  the  2 inch  plug  resulted  in  a doubling 
of  the  instrumentation  density. 

2.3.4.  Probe  Drive  Mechanism 

> 

To  carry  out  detailed  flow  field  measurements  in  the  region  of  interest, 
a probe  drive  system  was  built  and  is  attached  to  its  own  7 inch  plug  that  fitted 
[ into  the  eccentric,  rotatable  system.  A photograph  and  a line  drawing  of  the 

assembly  are  presented  as  Figs.  11  and  12,  respectively.  This  probing  system 
is  different  from  the  one  used  in  the  earlier  work  (Ref.  4),  although  some 
parts  are  the  same.  The  earlier  work  was  performed  with  a probe  system  that 
i extended  all  the  way  across  the  tunnel  to  measure  the  boundary  layer  interaction 

on  the  wall  opposite  to  the  one  on  which  the  probe  drive  was  mounted.  Since  the 
region  of  interest  extended  into  the  streamwise  corner  formed  by  the  test  wall 

i 

and  the  shock  generator,  the  earlier  system  would  be  expected  to  cause  local 
flow  blockage  if  the  main  streamlined  probe  support  were  mounted  close  to  the 
shock  generator.  Therefore,  a new  mechanism  was  designed  with  the  probe  ex- 
; tendina  from  the  wall  on  which  the  interaction  is  to  be  measured. 
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The  assembly  can  move  a probe  normal  to  the  wall  and  also  yaw  it  around 


its  tip  simultaneously.  This  was  made  possible  by  locating  the  probe  stem 
0.75  inches  off  center  in  a 2 inch  plug  that  is  free  to  rotate.  To  obtain  a 
good  bearing,  the  2 inch  plug  was  fabricated  from  stainless  steel  and  matched 
to  an  aluminum-bronze  bushing  in  the  larger  plug.  The  rotation  of  the  probe 
and  the  2 inch  plug  was  provided  by  a reversible  d.c.  motor  with  an  integral 
reduction  gear  box.  A worm,  mounted  on  a short  shaft  connected  to  the  motor 
by  a set  of  bevel  gears,  turned  the  worm  wheel  that  is  doweled  into  the  main 
probe  holder  shaft,  which  rotates  the  2 inch  plug.  This  probe  holder  shaft 
has  a cam  that  can  slide  up  and  down  in  a slot  to  allow  the  movement  of  the 
probe  formal  to  the  wall. 

The  translation  of  the  probe  was  achieved  by  moving  the  complete  yaw 
actuator  with  a precision  screw  drive.  This  screw  was  connected  by  a coupling 
to  a second  reversible  d.c.  motor. 

The  vertical  movement  was  measured  by  a precision  linear  transducer, 
which  was  calibrated  against  a dial  gage.  A ten  turn  potentiometer  connected 
by  a reduction  gear  to  the  small  shaft  of  the  worm,  indicated  the  angular  dis- 
placement. It  was  found  essential  to  be  able  to  repeat  very  accurately  the 
angular  settings  of  the  probe  with  respect  to  a fixed  reference,  even  after 
the  assembly  had  been  relocated  by  turning  it.  This  objective  was  achieved 
by  mounting  a small  mirror  on  top  of  the  probe  holder  shaft  and  using  the  re- 
flected image  of  a small  helium-neon  laser  beam  at  a distance  of  4 feet  as  an 
angle  indicator. 

Five  different  probes  that  can  be  mounted  in  the  probe  holder  sh art, 
were  constructed,  see  Fig.  13.  These  probes  were  sealed  with  an  0-ring  against 
the  2 inch  plug.  The  probes  will  be  discussed  in  the  next  three  sections. 
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2.3.5.  Cobra  Probes 


2 . 3 . 5 . 1 . Yaw  Angle  - Pitot  Pressure  Probe 
All  probes  were  fabricated  from  stainless  tubing  of  various  sizes. 

The  probe  stem,  made  of  two  concentric  tubes  silver  soldered  together,  has  a 
final  outside  and  inside  diameter  of  0.125  and  0.074  inches  respectively,  to 
provide  the  structural  stiffness  required  for  its  maximum  extension  of  2.1 
inches  from  the  wall.  These  stems  were  not  streamlined,  like  the  cobra  probe 
used  in  a similar  study  by  Peake  (Ref.  8),  since  large  flow  angle  variations 
with  distance  from  the  wall  are  present  and  can  be  expected  to  result  in  a 
greatly  enlarged  apparent  size  of  a streamlined  probe  stem  to  the  cross  flow, 
resulting  in  extraneous  flow  interference. 

The  yaw-pitot  probe  tip,  extending  0.75  inches  upstream  of  the  stem,  was 
manufactured  from  a set  of  three  hypodermic  tubes  soldered  together  for  rigidity. 
These  three  tubes  had  an  original  outside  diameter  of  0.024  inch,  but  were 
rolled  flat,  resulting  in  a tip  height  of  0.016  and  a width  of  0.080  inch.  The 
outer  tubes  were  cut  back  at  45  degrees  to  form  a cobra  probe,  of  which  the 
center  tube  functions  as  a rectangular  pitot.  A frontal  view  of  the  actual 
probe,  photographed  through  a microscope  after  all  the  tests  had  been  completed, 
is  shown  at  37  times  true  size  in  Fig.  14.  The  three  small  tubes  were  stepped 
up  to  an  outside  diameter  of  0.033  inches  and  led  through  the  probe  stem. 

Although  the  angle  sensitivity  of  this  yaw  probe  has  been  calibrated, 
as  will  be  discussed  in  Section  2.4.5.,  it  was  mainly  operated  as  a null  device. 
The  probe  was  moved  continuously  through  the  flow  at  a varying  extension  rate, 
while  being  kept  aligned  by  observing  the  differential  pressure  of  the  outer 
rubes.  This  resulted  in  a total  time  for  a 2 inch  survey  of  2 to  4 minutes, 
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depending  upon  the  skewness  of  the  profile  obtained.  This  extension  rate  was 
sufficiently  slow  so  that  no  pressure-time  lag  occurred. 

Periodic  checks  of  the  null  position  of  the  probe  were  made  by  locating 
it  in  the  free  stream  upstream  of  the  interaction  region  with  an  undeflected 
shock  generator.  A hysteresis  of  0.16  degrees  was  observed  in  the  null  posi- 
tion of  the  probe  if  this  angular  setting  was  approached  from  different  direc- 
tions. To  avoid  any  errors,  it  was  necessary  to  approach  the  null  from  the 
same  direction  in  any  particular  survey,  if  possible.  Yaw  angle  surveys  that 
have  a change  of  sign  in  the  angle  variation  with  distance  from  the  wall  have 
been  corrected  for  the  hyteresis  error  of  0.16  degrees  a posteriori. 

2 . 3 . 5 . 2 . Pitch  Angle-Pitot  Pressure  Probe 

To  investigate  the  boundary  layer  on  the  surface  of  the  shock  generator, 
it  was  found  important  to  fabricate  a cobra  probe  that  is  angle  sensitive  in  a 
plane  rotated  90°  from  the  one  of  the  yaw-pitot  probe.  The  drive  mechanism 
used  for  this  probe  is  the  same  as  for  the  yaw-pitot  probe.  The  operating 
procedure,  however,  had  to  be  different,  since  the  drive  cannot  vary  the  pitch 
angle  of  a probe,  thus  eliminating  the  possibility  of  using  this  probe  as  a 
nulling  device. 

The  movement  of  the  probe  tip  normal  to  the  surface  of  the  shock  gener- 
ator was  made  possible  by  mounting  the  probe  in  the  probe  holder  shaft  with  a 
180°  difference  in  orientation,  with  a probe  tip  location  1.5  inches  away  from 
the  center  of  rotation,  and  subsequently  yawing  the  mechanism.  The  pitch  angle 
of  the  flow  was  derived  by  using  the  angle  calibration  obtained  for  the  yaw- 
pitot  probe. 

The  minimum  outside  dimension  of  the  tip  is  0.011  inches  as  shown  in 
Fig.  15.  This  small  size  of  the  tip  resulted  in  a slower  time  response  to  a 
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pressure  change,  which  made  it  more  practical  to  measure  the  profiles  point 
by  point,  rather  than  through  continuous  movement.  An  enlarged  photograph  of 
the  tip  is  shown  in  Fig.  14. 

2.3.6.  Static  Pressure  Probe 

To  obtain  local  values  of  the  static  pressure  in  the  flow  field,  a small 
cone-cyl i rider  probe  was  designed,  to  be  used  in  conjunction  with  the  drive 
mechanism.  The  design  consists  of  a 20°  half-angle  cone  tip  on  a cylindrical 
body  constructed  from  0.025  inch  hypodermic  tubing.  Two  static  pressure  ori- 
fices of  0.008  inch  diameter,  180°  apart,  11  diameters  downstream  of  the  tip, 
and  followed  by  a 10  diameter  length  of  straight  tubing,  are  shown  in  Fig.  15. 
These  orifices  are  located  at  the  center  of  rotation  of  the  probe.  Photographs, 
taken  through  a microscope,  of  both  the  tip  and  the  static  pressure  orifices 
are  shown  in  Fig.  16. 

The  probe  was  traversed  continuously,  while  being  kept  aligned.  This 
procedure  involves  an  X-Y  plotter,  which  displays  the  position  and  the  angular 
orientation  of  the  probe.  The  probe  was  traversed  in  such  a way  that  the 
plotter  pen  followed  an  already  drawn  line  that  represented  the  angle  informa- 
tion obtained  from  the  yaw-pitot  probe.  This  procedure  is  straightforward  but 
tedious. 

Some  probe-wall  interference,  angle  sensitivity  and  probe  to  oblique 
shock  wave  interaction  has  been  observed  and  will  be  discussed  in  Section  2.4.6 
under  data  uncertainties. 

2.3.7.  Total  Temperature  Probes 

Two  fine  wire  probes  were  constructed,  according  to  a design  based  on 
hot-wire  concepts,  to  measure  the  local  total  temperature  through  the  inter- 
action flow  field.  The  probes  basically  consist  of  a fine  thermocouple  wire 
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spanning  two  supports.  These  probes  can  be  used  close  to  the  test  surface  by 
making  the  supports  at  a slight  angle  with  respect  to  that  surface  and  the 
wire  parallel  to  it.  Two  probes  of  this  type  were  made,  one  with  the  wire 
parallel  to  the  tunnel  wall,  indicated  as  total  temperature  probe  I in  Fig.  13, 
and  one  with  the  wire  parallel  to  the  shock  generator  surface,  total  temperature 
probe  II. 

The  actual  sensing  element  is  a chromel -al umel  thermocouple  junction  mid- 
way between  the  two  supports,  which  were  also  made  of  chromel  and  alumel.  All 
joints  were  welded  for  strength  purposes.  The  two  supports  pass  through  a 
ceramic  insulator  held  in  the  probe  stem  with  epoxy.  The  diameter  of  the  fine 
wire  is  0.002  inch  and  the  overall  aspect  ratio  was  50,  as  shown  in  Fig.  15. 

Two  photographs  of  the  tips  of  one  of  the  probes  are  presented  as  Fig.  17  at 
approximately  60  times  true  size.  One  can  observe  the  thermocouple  junction 
in  the  fine  wire  and  can  also  notice  that  the  support  tips  have  been  cut  back 
to  form  a wedge. 

The  total  temperature  surveys  were  performed  using  the  same  procedure  as 
was  used  for  the  static  pressure  probe  in  keeping  the  probe  aligned  with  the 
known  flow  direction  while  making  a traverse. 

2.3.8.  Data  Acquisition  Systems 

* 

During  the  course  of  the  study  many  different  configurations  of  the 
electronic  data  gathering  system  were  employed,  as  was  also  reported  in  the 
earlier  study  (Ref.  4).  The  basic  variation  is  the  rate  at  which  the  data 
was  taken,  the  system  involving  the  Scanivalve  being  the  slowest,  and  the  heat 
transfer  studies  having  the  highest  data  rate.  The  major  improvement  compared 
to  the  earlier  work,  concerned  the  speed  of  the  ITI  multiplexer,  which  has  been 
increased  from  100  to  200  readings  per  second  in  order  to  be  able  to  evaluate 
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the  tenperature -time  gradients  in  the  heat  transfer  technique  more  accurately. 
Periodic  checks  were  made  by  comparing  standard  inputs  with  the  final  output 
at  the  computer  end  to  insure  that  the  complete  system  operated  up  to  its 
specifications.  A typical  flow  diagram  showing  a system  configuration  for 
flow  field  surveys,  is  shown  as  Fig.  18  for  completeness. 

2 . 4 Estimated  Uncertainties  and  Reduction  of  Data 

2.4.1.  Transducer  Calibrations 

All  the  pressures  in  the  study  were  measured  by  Pace  variable  reluctance 
transducers  referenced  to  vacuum.  These  transducers  have  interchangeable  dia- 
phragms and  their  pressure  operating  range  can  thus  be  changed.  They  were 
calibrated  against  accurate  iieise  and  Wallace  & Tiernan  gauges. 

The  thermocouple  signals  were  fed  into  high  impedance  Datel  linear  ampli- 
fiers. The  calibration  was  made  by  switching  a standard  millivolt  source  to  the 
input  of  the  amplifiers  and  using  the  temperature -mi  11 i vol t tables  issued  by 
the  IJ.  S.  National  Bureau  of  Standards. 

The  traverse  distance  and  the  angular  displacement  of  the  probe  were 
indicated  by  a linear  transducer  and  ten  turn  potentiometer  respectively, 
which  were  calibrated  against  dial  indicators. 

All  the  transducers  were  calibrated  before  each  run.  Great  care  was 
taken  in  obtaining  precise  calibrations;  they  were  observed  to  be  linear  and 
repeatable  within  the  accuracies  quoted  in  Table  IV. 


TABLE  IV:  TRANSDUCER  ACCURACIES 


Transducer 

Full  Scale 

Accuracy 

% of  F.S . 

stagnation  pressure 

200  psia 

+0.4  psia 

+ 0.2 

pitot  pressure 

50  psia 

+0.1  psia 

+ 0.2 

surface  static  pressure 

15  psia 

+P .02  psia 

+ 0.14 

flow  static  pressure 

10  psia 

+0.04  psia 

+ 0.4 

Ap  pressure 

5 psia 

+0.005  psia 

+ 0.1 

temperatures 

530  °R 

+1.0  °R 

+ 0.2 

angular  displacement:  a 

o 

O 

<o 

+0.05° 

+ 0.08 

probe  traverse  distance:  Z 

2.1  inch 

_ _ 

+0.001  inch 

+ 0.05% 

2.4.2.  Shock  Generator  Deflection  Angles 

The  shock  generator  was  moved  manually  by  means  of  a screw  drive  manu- 
factured to  close  tolerances  (see  Fig.  4).  Its  deflection  angle  is  directly 
related  to  a mechanical  counter  geared  to  the  drive  shaft  and  the  various  de- 
flection angle  settings  could  be  repeated  to  within  0.05  degrees. 

Two  shock  generator  deflection  angles  are  defined  in  this  study.  One 
is  the  geometrical  generator  angle,  <j>g  , which  represents  the  physical  position 
of  the  shock  generator  surface  with  respect  to  the  tunnel  axis.  A second  is 
the  effective  shock  generator  angle,  aG  , which  is  a number  computed  from  the 
average  value  measured  by  the  static  pressure  taps  on  the  generator  and  the 
upstream  free  stream  Mach  number  using  oblique  shock  theory.  It  should  be 
noticed  here  that  these  effective  shock  generator  angles  will  slightly  depend 
upon  the  particular  location  of  the  generator  taps  since  it  can  be  expected 
that  the  generator  surface  static  pressure  distribution  is  not  uniform.  The 


"nominal  shock  generator  angle"  has  been  adopted  as  the  nomenclature,  indi- 
cating the  effective  or  geometrical  one  rounded  off  to  an  integer  value. 


The  displacement  effect  of  the  boundary  layer  growing  on  the  surface 
of  the  generator  is  considered  to  be  a first  order  approximation  of  the  dif- 
ference between  the  effective  and  geometrical  shock  generator  angles.  This 
displacement  thickness,  as  will  be  discussed  later,  was  measured  at  X = 7.6 
inches  for  nominal  generator  angles  of  4 and  10  degrees  and  was  found  to  be 
0.0287  and  0.0269  inches  respectively.  This  allows  one,  by  assuming  linear 
growth  of  the  displacement  thickness,  to  compute  (a^  - 4>g)  for  4°  and  10°  and 
to  compare  it  with  the  actual  difference  found.  The  results,  for  all  the  shock 
generator  angles  considered  in  this  study,  are  presented  in  Table  V. 


TABLE  V:  SHOCK  GENERATOR  ANGLES  (DEGREES) 


Nominal 
S.G.  Angle 

<$>G 

Geometrical 

aG 

Effecti ve 

(“G-^G) 

Measured 

(“G-4>g)  ^ 

Computed  from  6 

2 

1.65 

2.00 

0.35 

4 

3.75 

4.00 

0.25 

0.22 

6 

5.65 

6.10 

0.45 

8 

7.65 

8.00 

0.35 

10 

9.75 

10.00 

0.25 

0.20 

12 

11.80 

12.05 

0.25 

13 

12.80 

13.10 

0.30 

13.80 


14.10 


0.30 


2.4.3.  Surface  Static  Pressure 


The  most  usual  method  of  measuring  surface  static  pressure  (also 
employed  in  this  study)  by  using  a hole  drilled  in  and  normal  to  the  test 
surface,  is  subject  to  errors.  Great  care  has  been  taken  to  avoid  burrs, 
rounded  edges  or  other  imperfections,  but  even  when  the  hole  has  a near  per- 
fect geometry  a residual  error  remains.  Franklin  and  Wallace  (Ref.  9)  have 
measured  this  residual  "static-hole"  error  in  a most  definitive  study  by  using 
transducers  mounted  flush  with  the  surface. 

We  may  compute  this  "static-hole"  error,  using  the  universal  relation 
presented  in  Ref.  9 and  extending  it  to  compressible  flow  by  evaluating  the 
density  and  kinematic  viscosity  at  the  wall  temperature.  The  resulting 
"static-hole"  error  for  the  initial  boundary  layer  flow  with  a hole  diameter 
of  0.032  inch  is  found  to  be  10%  of  the  local  shear  stress.  But  since  the 
shear  stress  is  only  0.7%  of  the  local  static  pressure  level,  this  does  not 
produce  any  significant  error  even  if  the  shear  stress  is  allowed  to  increase 
by  a factor  of  five  in  the  interaction  area.  It  is  interesting  to  note  that 
this  result  does  not  directly  depend  upon  the  relative  size  of  the  hole  com- 
pared to  the  spatial  scale  of  the  boundary  layer. 

There  will  be  a second  error  in  the  static-hole  reading  if  pressure 
gradients  are  present,  which  the  Franklin  and  Wallace  correlation  did  not  take 
into  account.  This  error  can  also  be  expressed  in  terms  of  displacement,  e , 
of  the  effective  center  from  the  geometrical  center  and  can  accordingly  be 
estimated  by  making  the  physical  assumption  that  the  effective  center  remains 

1 G 1 

within  the  hole  diameter,  d (e.g.,  - ^ < j < >r)  • The  maximum  pressure  gra- 
dient in  the  present  work  has  been  observed  to  be  70%  per  inch  and  this  results, 
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with  a maximum  hole  diameter  of  0.032  inch,  in  an  upper  bound  of  the  pressure 
gradient  related  static-hole  error  of  + 1%. 

Additional  uncertainties  in  the  data  obtained,  arise  from  the  geometry 
of  the  test  surface.  Examples  of  these  would  be  leaks  along  the  edges  of  the 
2 inch  static  pressure  plug,  however  small,  which  could  not  be  eliminated,  and 
the  pressure-time  lag  in  the  Scani valve  (see  Section  2.3.2).  These  two  items 
can  be  appraised  if  one  realizes  that  any  one  static  pressure  trace  consists 
of  at  least  four  different  plug  locations,  each  sampled  with  opposite  directions 
in  Scanivalve  rotation. 

The  overall  uncertainty  of  the  static  pressure  data  is  estimated  to  be 
less  than  +2%  of  the  upstream  static  pressure  level. 

2.4.4.  Heat  Transfer  Studies 

2. 4. 4.1.  Heat  Transfer  Data  Reduction 

The  procedure  for  taking  the  actual  measurements  has  been  described 
in  Section  2.3.3  and  typical  temperature-time  traces  have  been  presented  as 
Fig.  10. 

Let  t indicate  the  time  at  which  the  slug  and  tunnel  wall  have  the 
o 

same  temperature  and  evaluate  the  heat  transfer  rate  at  time  tQ  . From  an 
exponential  curve,  fitted  to  a total  of  N points  symmetrically  distributed 
around  time  tQ  , one  can  evaluate  the  temperature-time  gradient  of  the  slug 
as  given  by  Equation  1. 
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where  i is  a successive  time-index  of  data  points  in  a trace. 


The  conservation  of  energy  requires  that  the  rate  decrease  of  the 
internal  energy  of  the  slug  is  equal  to  the  sum  of  the  convection,  radiation 
and  conduction.  By  neglecting  radiation  and  conduction  one  can  compute  the 
convective  heat  transfer  rate  directly  from  Equation  2: 


Vs  ,dTs. 

A$  (dt> 


Finally,  the  heat  transfer  coefficient  follows  from: 


(2) 


(3) 


It  should  be  noted  here  that  this  heat  transfer  coefficient,  as  it  is  evaluated 
within  the  interaction  region,  has  been  non-dimensional ized  by  the  mass  flow 
rate,  u-|  , in  the  free  stream  ahead  of  the  interaction.  This  definition 
also  incorporates  the  assumption  that  q_  is  proportional  to  p,  u,  (T  -T  ) 
for  small  test  to  test  variations  in  tunnel  total  temperature  and  wall  tempera- 
ture. This  means  that  if  P-|  » u-|  , Tw  and  Tgw  are  evaluated  at  the  specific 
conditions  of  each  individual  measurement  point,  Ch  will  be  invariant  for 
small  changes  in  test  conditions.  This  procedure  has  also  been  followed  in  the 
earlier  work  (Ref.  4) . 

Equation  (1)  has  been  evaluated  with  N = 18  in  this  study  and  since  the 
time  between  data  points  in  a trace  is  50  milliseconds,  the  exponential  curve 
spans  a time  period  of  0.85  seconds.  To  check  whether  the  results  are  consis- 
tent with  the  data  points  from  which  they  are  derived,  a second  data  reduction 
scheme,  replacing  Equation  (1),  has  also  been  employed.  This  second  method 
divided  the  set  of  N = 18  points  into  12  subsets  of  7 consecutive  points  each. 
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For  each  subset  the  temperature-time  derivative  was  evaluated  by  fitting  a 
straight  line  through  the  7 points.  The  heat  transfer  coefficient  was  obtained 
from  Equations  (2)  and  (3),  in  which  the  slug  temperature  at  the  central  point 
of  each  subset  of  7 was  substituted  for  the  wall  temperature.  The  final  heat 
transfer  coefficient  was  obtained  by  summing  the  12  coefficients  obtained  for 
the  subsets  and  dividing  by  12. 

This  second  procedure  was  used  to  screen  the  coefficients  derived  with 
the  use  of  Equation  (1)  for  consistency  with  the  original  data,  i.e.  if  the 
two  methods  produced  coefficients  that  were  not  within  5%  of  each  other,  the 
trace  was  discarded.  This  requirement,  however,  turned  out  to  be  satisfied 
for  99%  of  the  data,  proving  the  consistency  of  the  results  obtained. 

2 . 4 . 4 . 2 . Adiabatic  Wall  Temperature 

To  measure  the  local  adiabatic  wall  temperature  in  the  interactions 
one  needs  to  satisfy  some  very  stringent  requirements.  The  principal  one  is 
that  the  device  used  to  measure  it  should  not  change  the  conditions  of  heat 
transfer  that  would  exist  if  it  were  absent.  This  eliminates  all  conventional 
instrumentation  that  indicates  the  adiabatic  wall  temperature  directly.  An 
often  employed  alternative  is  to  measure  heat  transfer  rates  at  many  different 
driving  potentials  (Tw_Taw)  and  parameterizing  the  resulting  heat  transfer 
coefficients  with  respect  to  the  unknown  T . The  local  adiabatic  wall 
temperature  can  then  be  determined  from  the  requirement  that  the  heat  transfer 
coefficient  shall  be  invariant  for  various  (T  -T  ) . This  method,  however, 

W aW 

could  not  be  used  because  the  driving  potential  is  a characteristic  of  the 
present  facility  and  cannot  be  changed  at  will  for  a tunnel  wall  mounted  model. 
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A completely  different  and  simple  way  of  solving  this  potential  problem 
is  to  discard  the  idea  of  a local  adiabatic  wall  temperature  in  the  inter- 
action entirely  as  done  in  this  study. 

Measurements  of  recovery  factors,  r , have  been  obtained  for  supersonic 
flows  on  insulated  surfaces  with  constant  pressure  and  have  been  reported  to 
be  approximately  0.89  where  r is  defined  in  Equation  (4): 


r 


(4) 


This  can  be  used  to  write  the  adiabatic  wall  temperature  directly  in  terms  of 
stagnation  pressure  for  a fixed  Mach  number: 

^rMl-Df  (5) 

o o 

Tco 

Hence,  with  r = 0.89  , M = 2.94  and  = 0.3665  this  results  in 

o 

Ta„  = 0.930  T0  (6) 


Equation  (6)  is  valid  for  the  initial  constant  pressure  turbulent  boun- 
dary layer  in  this  study.  This  relation  for  T&w  , evaluated  with  TQ  for 
every  specific  test,  has  been  u.:ed  throughout  this  study.  If  used  in  Equation 
(3)  it  results  in  a concise  definition  of  the  heat  transfer  coefficient,  since 
the  coefficient  becomes  directly  proportional  to  the  dimensional  heat  transfer 
rate  in  nonadiabatic,  physical  situation  and  no  computations  involving  local 
adiabatic  wall  temperatures  have  to  be  executed. 
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2.4.4. 3.  Van  Priest  I!  Skin  Friction  Formula 
It  is  recommended  by  Hopkins  and  Keener  (Ref.  1)  that  the  van  Driest  II 
correlation  with  a Reynolds  analogy  factor  of  1.2  be  used  to  calculate  the 
heat  transfer  coefficient  for  near  adiabatic  tunnel  walls  exposed  to  super- 
sonic Mach  numbers.  This  scheme  is  claimed  to  predict  the  heat  transfer  coef- 
ficient within  + 10%  for  near  adiabatic  temperature  conditions.  A summary  of 
the  scheme,  presented  as  a transformation  to  an  incompressible  situation  (Ref.  1), 
is  given  in  the  following  equations. 


Ree  ' ^ Ree 


1/Cf  = 17.08  (log-|Q  Re0)  + 25.11  1 og-j ^ Re^  + 6.012 


Karman-Schoenherr  Eq. 


Cf  = VFc 


F„  = ("9—  r Mj7(sin-1  a'  + sin"1  b') 


a'  = (2A2-B)/(4A2+B2) 


B’  = B/(4A2+B2)1/2 

a2  = 
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B = (HrY-C-T*) 


Chv0  = ? Cf  = 0‘6  Cf 


The  heat  transfer  coefficients  in  the  boundary  layer  upstream  of  the  inter- 
action are  computed  for  the  range  of  test  conditions  specified  in  Table  III 
of  Section  2.2.3  and  using  the  following  data: 

pQ  = 100  (psia) 

0 = 0.025  inch 

r = 0.89 
s = 1.2 
M = M,  = 2.94 

oo  1 

T = T,  = 0.3665  T 

00  1 0 

= Vi-, 

resulting  in:  Ch  = 0.000711  +0.000012 

These  results  show  that,  theoretically,  the  variations  of  temperature  condi- 
tions amount  only  to  a + 1.7%  variation  of  the  heat  transfer  coefficient  in 
the  upstream  boundary  layer  flow. 

2. 4. 4. 4.  Heat  Transfer  Variation  with  Reynolds  Number 
Prior  to  making  the  heat  transfer  measurements  in  the  interaction  region 
at  100  psia  stagnation  pressure,  some  preliminary  measurements  were  carried  out 
in  the  boundary  layer  upstream  of  the  interaction  region  with  an  undeflected 
shock  generator  and  a variable  stagnation  pressure  ranging  from  75  to  200  psia. 
The  resulting  heat  transfer  coefficients  of  these  measurements  divided  by  the 
van  Driest  value,  computed  as  set  out  in  Section  2.4.4,  are  presented  in  Fig.  19 
as  a function  of  both  Reynolds  number  based  on  momentum  defect  thickness  as  well 
as  stagnation  pressure  in  psia.  One  observes  that  the  measurement  points  are. 


on  the  average,  10“  higher  than  the  predicted  value,  and  fall  within  a band 
between  -15%  and  +10%.  If  one  realizes  that  the  dimensional  heat  transfer 
rate  varied  from  0.6  to  2.4  (Btu/ft  sec.),  although  the  heat  transfer  coeffi- 
cient does  not  change  significantly  over  the  Reynolds  number  range,  and  that 
the  same  range  in  dimensional  heating  rates  has  been  observed  in  the  interaction, 
one  can  conclude  that  the  technique  employed  certainly  should  indicate  the  cor- 
rect trends  and  levels  of  in  the  interaction  region.  In  other  words,  one 
can  consider  the  correct  verification  of  the  change  in  with  Re0  to  be  a 
calibration  of  the  technique,  although  it  is  not  used  as  such. 

2 . 4 . 4 . 5 . Heat  Transfer  Data  Uncertainty 
Assumptions  of  the  technique,  e.g.  neglection  of  radiation  and  conduc- 
tion losses,  have  been  verified  and  found  to  be  correct.  A good  estimate  of 
these  errors  has  been  given  by  Reinecke  (Ref.  10)  and  will  not  be  repeated 
here.  The  largest  source  of  error  in  the  present  results,  although  improved 
compared  to  the  earlier  work  (Ref.  4),  is  found  to  be  in  the  basic  inaccuracies 
of  the  temperature  measurements,  (+1.0°R)  followed  by  a differentiation  of  the 
data.  The  final  heat  transfer  data  in  the  interaction  is  considered  to  have 
an  accuracy  of  +15%  of  the  van  Driest  value.  Averaged  values  of  the  results 
obtained,  however,  are  believed  to  be  much  better  than  the  number  quoted  above. 
2.4.5.  Cobra  Probe  Surveys 

2 . 4 . 5 . 1 . Cobra  Probe  Displacement  Effect 
A pitot  probe  or  a cobra  probe  placed  in  a shear  flow  does  not  measure 
the  true  pressure  that  exists  at  the  geometric  center  of  the  probe  in  a flow 
with  no  shear.  The  effective  center  is  said  to  be  "displaced"  from  the  geo- 
metric center.  If  the  magnitude  and  the  direction  of  this  displacement,  A , are 
known,  corrections  to  the  pitot-probe  measurements  can  be  made  for  this  effect. 
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Experiments  to  appraise  this  effect  have  been  made  by  Young  and  Maas 
(Ref.  11)  at  low  subsonic  speeds.  They  found  that  the  displacement  towards 
the  region  of  higher  velocity,  A , divided  by  tne  probe  height,  H , had  a con- 
stant value  of  0.18.  This  value  of  A/H  = 0.18  was  revealed  to  be  independent 
of  H , the  mean  pitot  pressure,  and  the  magnitude  of  the  transverse  total  head 
gradient  across  the  probe  face.  More  recently  Allen  (Ref.  12)  studied  this 
displacement  effect  in  a turbulent  boundary  layer  at  a Mach  number  of  2 and  his 
experiment  yielded  an  average  value  of  0.38  for  A/H  ; it  is  interesting  to 
note  that  Allen  also  observed  A/H  to  be  independent  of  H , pitot  pressure, 
and  its  gradient  across  the  probe  face.  An  experimental  study  by  Johannesen 
and  Maier  (Ref.  13)  at  a Mach  number  of  1.98,  however,  showed  that  there  was 
no  displacement  effect  in  supersonic  shear  flow. 

In  view  of  the  contradictory  findings  at  Mach  2 and  the  fact  that  the 
height  of  the  present  yaw-angle/pitot-pressure  probe  is  only  3.2%  of  the  initial 
boundary  layer  thickness,  no  correction  to  yaw-pitot  measurements  have  been 
made.  This  will  not  affect  the  uncertainty  of  the  measurements  in  general,  but 
may  be  important  for  the  measurement  points  within  one  probe  height  from  the 
test  wal 1 . 

2. 4. 5. 2.  Cobra  Probe  Angle  Sensitivity 
Although  the  yaw  angle-pitot  pressure  probe  has  mainly  been  used  as  a 
nulling  device,  it  was  nevertheless  calibrated  as  a direct  reading  device. 

This  calibration  was  achieved  by  moving  the  probe  off  the  null  position  and 
observing  the  difference  in  pressure  in  the  outer  tubes,  Ap  . The  procedure, 
performed  at  various  static  pressure  levels  and  pitot  pressures,  showed  that 
Ap  varied  linearly  with  the  probe  deflection  angle,  ap  , for  small  deflec- 
tions (-10°  < a < 10°).  The  data  obtained  can  be  correlated  according  to 
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Equation  (11)  as  shown  in  Fig.  20. 

where  C-j  = 0.036  + 0.004  (deg~^)  for  -10°  < Up  < 10°  and  1.1  < M < 3.0  . 
This  equation  simply  says  that  the  pressure  difference,  Ap  , is  proportional 
to  a and  is  also  proportional  to  the  difference  between  pitot  pressure  and 
static  pressure.  For  supersonic  flow  the  angle  sensitivity  can  also  be  ex- 
pressed with  Mach  number  dependence  only  by  using  the  Rayleigh  pitot  formula 
and  y = 1.4  . 

3.5  2.5 

(—1—)  -1)  (12) 
7M  -1 

1.1  < M < 3.0 

The  constant  C-|  is  primarily  a function  of  the  angle  at  which  the  outer  tubes 
are  cut  back  (45°  in  this  study)  and  the  ratio  of  the  probe  tip  width,  W , over 
the  height,  H . (^/h  = 6 for  the  yaw-pitot  probe.) 

2. 4. 5. 3.  Cobra  Probe  Data  Uncertainties 
As  discussed  in  Section  2. 3. 5.1.  the  yaw-pitot  probe  has  mainly  been 
operated  as  a null  device.  Small  misalignment,  however,  could  still  occur  and 
this  has  been  corrected  for  by  using  Equation  (11).  The  minimum  accuracy  of 
the  flow  angles  obtained  can  be  assessed  by  combining  the  accuracy  of  the  Ap 
transducer  with  the  angle  sensitivity  of  the  probe.  From  Equation  (12)  one 
obtains  for 

M = 1.1:  Ap/Aoip  = 0.041  p]  (deg-1) 

and  M = 3.0:  Ap/Acip  = 0.40  p-j  (deg-^) 
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The  minimum  static  pressure  is  3 psia;  this  combined  with  the  lowest  angle 
sensitivity  at  M = 1.1  gives: 

(Ap/Aa  ) = 0.041  * 3 = 0.123 

M mi  n y 

r 

With  an  accuracy  of  +0.005  (psia)  for  the  Ap  transducer,  there  results  a 
maximum  error  for  the  flow  angle  due  to  probe  misalignment  of  +0.04°. 

These  estimates  combined  with  the  a transducer  accuracy  of  + 0.05° 
and  the  repeatability  of  angular  settings  of  the  probe  drive  mechanism  results 
in  the  overall  estimate  that  the  flow  angles  measured  by  the  yaw-pitot  probe 
are  accurate  within  +0.1  (degrees). 

There  is  usually  no  problem,  as  there  is  with  a static  pressure  probe, 
in  obtaining  accurate  pitot  tube  measurements  across  an  oblique  shock  wave, 
since  the  pitot  pressure  is  influenced  by  the  bow  shock/oblique  shock  inter- 
action at  distances  less  than  one  tip  height  as  shown  by  Bannink  and  Nebbeling 
(Ref.  14).  The  overall  accuracy  of  the  pitot  measurements  is  therefore  esti- 
mated to  be  well  within  +1%  of  the  maximum  pitot  pressure. 

2.4.6  Static  Pressure  Surveys 

The  particular  types  of  static  pressure  probe  employed  in  this  study, 
the  cone-cylinder  type,  is  known  to  create  extraneous  interference  if  the 
probe  extends  across  an  oblique  shock  wave.  Back  and  Cuffel  (Ref.  15)  have 
appraised  this  interference  effect  between  the  probe  bow  shock  and  the  oblique 
shock  wave  at  a Mach  number  of  3.5  and  with  an  oblique  shock  wave  created  by 
an  8.3°  flow  deflection.  Errors  up  to  50%  of  the  static  pressure  difference 
across  the  shock  were  measured.  It  should  be  remembered,  however,  that  these 
errors  only  occur  if  the  probe,  which  is  approximately  20  diameters  long. 


extends  across  the  oblique  shock  wave  and  thus  only  affects  the  measurements 
in  the  vicinity  of  the  shock. 

A superior  method  of  making  static  pressure  measurements  in  the  vicinity 
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of  oblique  shock  waves  is  the  conical  static  pressure  probe  technique,  which 
employs  a static  hole  on  the  conical  surface  of  a probe  tip.  This  method  has 
the  advantage  that  the  static  pressure  can  be  measured  up  to  one  cone-tip 
length  from  the  shock,  similar  to  the  behavior  of  a pitot  probe  in  the  vicinity 
of  a shock.  This  conical  static  pressure  probe  technique  has  some  drawbacks, 
namely  an  extensive  Mach  number  calibration  is  needed,  and  the  probe  does  not 
read  the  static  pressure  directly.  An  iterative  procedure  is  then  required 
because  the  Mach  number  is  not  known  independently  from  the  static  pressure  to 
be  measured.  The  angle  sensitivity  of  this  probe,  as  reported  by  Norris  (Refs. 
16  & 17),  can  be  viewed  as  an  advantage  since  it  could  be  used  as  a yaw  meter, 
although  this  sensitivity  is  not  as  large  as  the  one  for  the  cobra  probe  em- 
ployed in  this  study.  This  conical  probe  has  not  been  employed  in  the  present 
study  because  of  its  limited  usefulness  compared  to  the  effort  required. 

The  angle  sensitivity  of  the  cone-cylinder  probe  has  been  determined  in 
the  constant  pressure  boundary  layer  upstream  of  the  interaction,  as  this  might 
give  an  insight  to  the  errors  caused  by  an  misalignment  that  might  occur  in  the 
operational  procedure  described  in  Section  2.4.6.  The  results  are  shown  in 
Fig.  21.  One  can  observe  that,  with  the  probe  aligned,  the  direct  reading  is 
within  +0.4%  of  the  correct  value,  which  compares  favorably  with  the  results 
obtained  by  Cronvich  (Ref.  18).  With  a misalignment  angle  within  +4°  the  probe 
still  records  the  correct  value  within  +1.6%.  Above  12°,  however,  the  reading 
rapidly  deteriorates. 
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The  cone-cylinder  probe  also  causes  a wall  interference  at  some  distance 
from  the  test  wall  as  the  tip  bow  shock  reflects  from  the  wall  and  impinges  on 
the  cylinder  in  the  neighborhood  of  the  static-holes.  This  problem  does  not 
arise  when  the  probe  body  is  sufficiently  close  to  the  wall  as  shown  in  Fig.  21, 
where  the  data  for  a static-hole  to  wall  distance  of  0.050  inch  is  shown;  in 
that  case  the  probe  reads  the  correct  pressure  as  indicated  by  a wall  surface 
static-hole. 

Other  uncertainties  in  the  static  probe  readings  remain  because  of  its 
long  body  (20  diameters),  and  thus  its  unknown  capability  to  resolve  continuous 
static  pressure  gradients,  and  its  unknown  response  to  strong  streamline  curva- 
ture. With  the  exclusion  of  the  data  obtained  within  the  vicinity  of  the 
oblique  shock  wave  and  the  data  that  shows  obvious  wall  interference,  the 
readings  of  the  cone-cylinder  static  pressure  probe  are  nevertheless  estimated 
to  have  an  accuracy  within  +5%  of  the  upstream  static  pressure,  p-j  . 

2.4.7.  Total  Temperature  Surveys 

The  fine  wire  probes  used  in  this  study  are  designed  according  to  hot- 
wire concepts  and  do  not  record  the  total  temperature  directly,  but  rather  indi- 
cate a wire  recovery  temperature,  which  is  equal  to  the  adiabatic  wire  recovery 
temperature  if  conduction  and  radiation  losses  can  be  neglected.  Even  if  these 
radiation  and  conduction  losses  cannot  be  neglected,  as  is  the  case  in  high 
temperature,  low  density  flows,  it  has  been  shown  by  Vas  (Ref.  19)  that  the 
local  total  temperature  can  still  be  determined  with  sufficient  accuracy  by 
such  a fine  wire  probe. 

The  adiabatic  wire  recovery  temperature,  as  reported  by  Dewey  (Ref.  20) 
is  equal  to  C.950  times  the  stagnation  temperature  at  Mach  3.  The  present 
experiments  have  shown  a wire  recovery  temperature  of  0.948  Tq  in  the  free 
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stream  and  thus  indicates  a loss  by  conduction  of  0.2%.  The  present  wire 
temperatures  have  been  reduced  to  local  total  temperatures  by  assuming  a con- 
stant recovery  factor  for  the  fine  wire.  The  recovery  factor  is  defined  as 
a nondimensional  enthalpy  increase: 


r 


Twi'T 


03) 


with  M = 2.94  , = 0.3665  , Twi  = 0.948  TQ  and  Tq  = T^  in  free  stream, 

o 

hence 


r = (V— )/(l-  X-)  = 0.918 
o o 

From  this  result  we  can  immediately  compute  the  total  temperature  as  a function 
of  Mach  number  from: 


X-  = (1  +0.2  M2)  with  y=1.4 
't 

and  from  Equation  (13) 


Substitute  (14)  in  (15) 


VTo 


r+(l-r)(l+0.2Nr) 


-1 


r = 0.918 


(14) 


(15) 


(16) 
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To  : tunnel  stagnation  temperature 

Tt  : local  total  temperature 

T • : wire  recovery  temperature 

M : local  Mach  number 

This  Mach  number  dependence  does  not  result  in  any  difficulties,  since  it  is 
already  known  from  the  pressure  measurements. 

It  has  been  assumed  that  the  variations  in  the  ratio  of  tunnel  wall 
temperature  to  tunnel  stagnation  temperature  do  not  influence  the  values  ob- 
tained for  T./T  in  the  interaction  flow  field.  This  is  based  on  the  obser- 

t o j 

vations  that  the  ratio  is  nearly  a constant  of  1.21  +0.04  . The  small 

aw 

variations  in  this  parameter  are  mainly  limited  to  the  laminar  sublayer,  which 
is  of  such  a small  physical  extent  from  the  wall  that  no  data  points  have  been 
obtained  in  it. 

The  total  uncertainty  of  the  total  temperature  results  is  estimated  to 
be  within  +1%  of  the  tunnel  stagnation  temperature. 

2.4.8.  Derived  Flow  Field  Variables 

With  the  measured  values  of  pitot  pressure,  static  pressure,  yaw  angle 
and  total  temperature  one  can  compute  the  Mach  numbers,  u and  v velocity 
components  and  the  densities  throughout  the  flow  field.  But  since  the  various 
measured  surveys  do  not  have  an  identical  number  of  points  in  the  Z-direction, 
one  has  to  set  up  a standard  array  of  Z values  at  which  the  measured  quan- 
tities are  to  be  known.  This  process  does  not  involve  any  data  interpolation 
since  continuous  surveys  have  been  taken  with  many  more  data  points  than  needed 
for  the  actual  standard  Z array. 

The  values  of  Mach  number,  u and  v velocity  components  and  density 


follow  from  the  following  equations: 
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Y^T 


~ = (1  + ^ M2)  if  ^ < 1 .893 
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if  jp>  1.893 


T = Tt/(1+  M2) 


p = P/RT 


u = MART  cos  a 


v = MART  sin  a 


(18) 

(19) 

(20) 
(21) 
(22) 


In  addition,  the  velocity  component  in  the  Z-di recti  on  can  be  deduced 
from  the  u,  v components  and  the  density  p by  integration  of  the  continuity 
equation : 


(fx  (pu)  + |y  (pv))  dz1 


(23) 


It  should  be  noted  that  the  actual  values  of  w obtained  should  be  small  com- 
pared to  ^ since  the  total  Mach  number  is  based  upon  the  u and  v 
velocities  only  and  could  be  considerably  in  error  if  large  relative  values 
of  w are  found.  This  requirement  can  be  checked  a posteriori  however.  Two 
pitch  angles  can  now  be  defined,  8 and  Y : 
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8 = tan"1  (|j) 
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tan 


( 


-w  ) 
/ u2+v2 ' 


(25) 


2.5  Uniformity  of  Initial  Boundary  Layer  Flow 

Most  of  the  previous  experiments,  Stanbrook  (Ref.  21),  McCabe  (Ref.  22), 
Lowrie  (Ref.  23),  Token  (Ref.  24)  and  Peake  (Ref.  25),  have  been  conducted  with 
the  turbulent  boundary  layer  developing  on  the  sidewall  of  a two-dimensional 
nozzle.  This  sidewall  boundary  layer  will  be  transversely  non-uniform  due  to 
the  presence  of  transverse  pressure  gradients  in  the  expanding  part  of  the  noz- 
zle which  causes  the  boundary  layer  flow  in  the  nozzle  to  accumulate  near  the 
center  of  the  sidewall.  The  magnitude  of  this  undesired  effect  increases  as 
the  nozzle  Mach  number  increases. 

In  some  of  these  studies  the  non-uniformity  has  been  reduced  by  taking 
appropriate  counter  measures.  Lowrie  re-designed  the  nozzle  layout.  Peake 
installed  a V-type  vortex  generator  at  Mach  4 to  "smear  out"  the  non -uniformity, 
but  was  only  moderately  successful  because  large  transverse  non-uniformities 
still  remained  which  at  Mach  4 reportedly  are  of  the  order  of  +8%  in  static 
pressure  and  +50%  in  displacement  and  momentum  defect  thickness.  The  best 
counter  measure  has  probably  been  taken  by  Token  who  only  used  the  central  33% 
of  the  sidewall  as  a test  region  but  still  was  left  with  a 13  inch  wide  area. 

In  that  study,  however,  transverse  pressure  differences  of  +3.5%  remained  in 
the  test  region.  These  non-uniformities  basically  limit  the  experimental  accu- 
racies that  can  be  obtained  within  the  shock  wave/boundary  layer  interaction. 

Recognizing  this  potential  problem,  the  test  section  and  shock  gener- 
ator have  been  designed  such  that  the  boundary  layer  entering  the  interaction 
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originates  from  the  contoured  wall  in  the  two-dimensional  nozzle,  as  is  shown 
in  Fig.  22.  This  boundary  layer  has  been  surveyed  extensively  and  shown  to  be 
extremely  uniform  transversely  across  the  central  6 inches  of  the  8 inch  wide 
test  section.  The  magnitude  of  the  uniformity  is  measured  to  be  +0.7%  in 
static  pressure  and  + 5%  in  boundary  layer  displacement  and  momentum  defect 
thicknesses  (Fig.  23).  This  level  of  uniformity  is  considered  to  be  an  impor- 
tant improvement  over  previous  studies  in  which  flow  field  yaw  angles  have 
been  measured. 


2.6  Summary  of  Estimated  Uncertainties  of  Data 

The  final  experimental  uncertainties,  as  discussed  in  Section  2.4,  are 
larger  than  the  accuracies  of  the  transducers  used  to  obtain  the  various  data. 
This,  however,  does  not  imply  that  the  uncertainties  are  inherent  to  the  basic 
instrumentation  techniques  employed.  The  specific  application  of  the  tech- 
niques in  each  part  of  the  interaction  is  the  factor  which  determines  the 
final  uncertainty  of  the  data.  The  percentages  given  are  therefore  an  upper 
bound  for  the  accuracy  and  many  results  are  believed  to  be  more  accurate  than 
the  numbers  quoted  in  Table  VI. 


TABLE  VI:  SUMMARY  OF  DATA  UNCERTAINTIES 


Quantity  Measured 

Symbol 

Uncertainty 

stagnation  pressure 

Po 

+0.5%  at  100  psia 

stagnation  temperature 

To 

+0.5%  of  450°R 

surface  static  pressure 

P 

+2%  of  Pi 

surface  heat  transfer 

Ch 

+1 5%  of  Chi 

flow  static  pressure 

P 

+5%  of  Pi 

pitot  pressure 

Pt 

+1%  of  p^i 

local  total  temperature 

Tt 

+1%  of  T0 

flow  yaw  angle 

a 

+0.1  degree 

probe  location 

Z 

+0.002  inch 

effective  shock 

generator  angle 
. ■ - . ■ . — ■ ■ — . — — - 

aG 

+0.1  degree 

2.7  Summary  of  Data 

The  data  obtained  is  graphically  presented  in  10  detailed  appendices, 

A through  J. 

Appendix  A contains  the  surface  static  pressure  distributions  obtained 
along  instrumentation  lines  parallel  to  the  Y axis  for  various  shock  generator 
angles  and  X locations.  Eoth  the  computed  shock  location  as  well  as  the  static 
pressure  level  measured  at  the  surface  of  the  shock  generator  are  indicated  for 
each  trace. 

Appendix  B shows  the  isobar  patterns  obtained  directly  from  the  data  in 
Appendix  A.  The  pressure  p-|  and  p£  are  the  static  pressures  in  front  of 

I 

the  interaction  and  on  the  surface  of  the  generator.  The  increase  in  surface 
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static  pressure,  going  from  p^  to  p2  , is  divided  into  10  equal  steps  and 

P-Pi 

the  corresponding  10  isobars  at  constant  values  of  (- — ' ) are  shown  in 

p2"pl 

the  figures  of  the  appendix. 

Appendix  C summarizes  the  heat  transfer  results,  non-dimensional ized 
by  the  computed  van  Driest  II  values,  obtained  for  the  initial  boundary  layer 
(see  Section  2. 4. 4. 3).  The  data  shown  is  taken  at  instrumentation  stations 
parallel  to  the  Y axis  at  various  X locations  and  for  7 shock  generator  angles 
ranging  from  2°  to  14°  in  2°  increments.  The  corresponding  surface  heat 
transfer  patterns  are  presented  in  Appendix  D. 

Appendices  E and  F contain  the  results  obtained  with  the  yaw-angle/ 
pitot-pressure  probe.  The  tick  marks  on  the  Z = 2.0  inches  axis  indicate  the 
level  which  the  traces  theoretically  should  approach  at  large  distances  from 
the  wall.  The  Appendices  total  49  surveys  each,  which  can  be  divided  into  7 
groups  as  shown  in  Table  VII. 


TABLE  VII:  YAW-PITOT  SURVEY  GROUPINGS 


NSGA 

degrees 

X 

inches 

Number  of 
surveys 

0 

3.561 

5 

4 

7.579 

16 

4 

7.079 

5 

4 

5.079 

3 

10 

7.599 

14 

10 

7.099 

3 

10 


5.099 


3 


Appendix  G presents,  in  graphical  form,  the  results  of  the  flow  field 
measurements . The  data  affected  by  probe-wall  interference  (Section  2.4.6) 


has  been  eliminated  and  replaced  by  a linear  interpolation  of  the  values  mea- 
sured at  Z = 0.25  inch  and  by  the  surface  instrumentation.  Figure  G-8  shows 
two  profiles  where  the  cone-cylinder  static  pressure  probe  interfered  with 
the  shock-wave  at  larger  values  of  Z resulting  in  off  scale  values  which  have 
been  replaced  by  a constant  estimated  level.  A total  of  29  surveys  are  plotted 
which  can  be  divided  into  groups  as  shown  in  Table  VIII: 


TABLE  VIII:  STATIC  PRESSURE  SURVEY  GROUPINGS 


NSGA 

degrees 

X 

inches 

Number  of 
surveys 

4 

7.579 

12 

4 

5.079 

3 

10 

7.599 

11 

10 

5.099 

3 

The  data  obtained  with  the  fine-wire  thermocouple  probe  I is  plotted 
in  Appendix  H.  This  data  can  be  divided  into  similar  groups  as  the  data  of 
Appendix  G,  as  shown  in  Table  IX. 

The  final  two  Appendices,  I and  J,  contain  the  results  obtained  with 
the  pitch-pitot  probe  and  thermocouple  probe  II  respectively.  The  final  three 
traces  of  Appendix  J,  taken  at  the  same  location,  with  the  same  shock  generator 
angle,  but  with  different  shock  generator  temperatures  (see  Ref.  5),  show  the 
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limited  extent  to  which  the  shock  generator  temperature  influences  the  total 
temperature  distribution  in  the  interaction  flow  field. 


TABLE  IX:  TOTAL  TEMPERATURE  SURVEY  GROUPINGS 


m 

Number  of 
surveys 

0 

3.561 

1 

4 

7.579 

10 

4 

5.079 

3 

10 

7.599 

9 

10 

5.099 

3 

3.  CONCLUDING  REMARKS 


A detailed  study  of  three-dimensional  flow  fields  and  surface  pheno- 
mena caused  by  the  interaction  of  an  oblique  shock  wave,  a shock  generator 
and  a high  Reynolds  number  turbulent  boundary  layer,  was  carried  out.  This 
initial  turbulent  layer  has  been  shown  to  be  uniform  in  the  transverse  direc- 
tion across  the  central  part  of  the  test  wall.  The  main  results  obtained 
during  the  present  contract  effort  include  detailed  surface  static  pressure 
distributions  and  surface  heat  transfer  patterns  from  which  a physical  insight 
into  the  spatial  scale  of  the  phenomena  may  be  derived.  In  addition,  detailed 
flow  field  measurements  were  conducted  to  fully  resolve  the  flow  field,  in- 
cluding all  three  velocity  components.  This  combination  of  detailed  surface 
data  and  associated  flow  field  surveys  provides  the  opportunity  to  check  how 
the  two  sets  are  related  and  may  lead  to  a deeper  understanding  of  the  problem 
under  consideration. 

Special  attention  was  given  to  the  mathematical  formulation  of  the  physical 
boundary  conditions  which  existed  in  the  experimental  configuration.  In  addi- 
tion, the  validity  of  the  experimental  methods  employed  to  obtain  the  various 
data  was  discussed  extensively  and  was  complemented  by  a detailed  tabulation 
of  the  estimated  upper  bounds  for  the  accuracies  of  the  various  data.  These 
considerations  are  considered  to  be  important  since  the  final  analysis  of  the 
experimental  data  obtained,  will  heavily  rely  on  the  accuracy  achieved.  No 
discussion  in  terms  of  the  physical  interpretation  of  the  presented  data  in 
relation  to  the  problems  of  viscous-inviscid  interactions  is  given. 
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Figure  1.  Schematic  of  Experimental  Configuration 
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Figure  2.  Coordinate  System. 
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Figure  9.  Drawing  of  Slug  Calorimeter  Arrangemen 
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Figure  10.  Temperature-Time  History  of  Slug  Calorimeter. 
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Fioure  12.  Drawinc  of  Probp  Drive  Mechanism. 
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Figure  13.  Photograph  of  Five  Flow  Field  Probes. 
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Figure  17.  Photograph  of  Total  Temperature  Proba 
Supports  and  Thermocouple  Junction. 
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Figure  19.  Heat  Transfer  Versus  Reynolds  Number. 
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Figure  A-2:  Surface  Static  Pressure,  X = 5.1  inches. 
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Figure  A-4:  Surface  Static  Pressure,  X = 7.1  inches. 
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Figure  A-5:  Surface  Static  Pressure,  X = 7.6  inches. 
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Fiaure  B-2:  Surface  Isobar  Pattern,  4 decrees 


B-3:  Surface  Isobar  Pattern,  6 decrees. 


Figure  B-4:  Surface  Isobar  Pattern,  8 degrees. 


Figure  B-5:  Surface  Isobar  Pattern, 10  degrees 


CH/CHVD.  HEAT  TRANSFER  RATIO 


HOI* INAL  SHOCK 
GENERATOR  ANGLE 


DISTANCE  YG  (INCHES) 

Fiaure  C-3:  Surface  Heat  Transfer,  X = 7.1  inches. 


Surface  Heat  Transfer  Pattern,  4 decrees. 


Surface  Heat  Transfer  Pattern,  6 degrees. 


Figure  D-3:  Surface  Heat  Transfer  Pattern,  8 degrees 


Surface  Heat  Transfer  Pattern,  10  denrees. 
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Figure  D-5:  Surface  Feat  Transfer  Pattern,  12  degrees. 
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Figure  E-l:  Yaw  Angle  Surveys,  X = 3.56  inches 
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Figure  E-3:  Yaw  Angle  Surveys,  X = 7.6  inches 
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Figure  E-6:  Yaw  Angle  Surveys,  X = 7.6  inches. 
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Figure  E-8:  Yaw  Angle  Surveys, 
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Figure  E-10:  Yaw  Angle  Surveys,  X = 7.6  inches 


Figure  E-ll:  Yaw  Angle  Surveys,  X = 7.6  inches. 
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Figure  E-14:  Yaw  Angle  Surveys,  X = 7.1  inches 
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Figure  E-15:  Yaw  Angle  Surveys,  X = 5.1  inches. 
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Fiaure  F-3:  Pitot  Pressure  Surveys,  X = 7.6  inches 
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Figure  F-4:  Pitot  Pressure  Surveys, 
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Fiaure  F-7:  Pitot  Pressure  Surveys,  X = 7.1  inches 


Figure  F-9:  Pitot  Pressure  Surveys,  X = 5.1  inches. 
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Figure  F-11:  Pitot  Pressure  Surveys,  X = 7.6  inches. 
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Figure  F-15:  Pitot  Pressure  Surveys,  X = 5.1  inches, 


Static  Pressure  Surveys,  X = 7.6  inches, 


Figure  G-3:  Static  Pressure  Surveys,  X = 7.6  inches. 


Pressure  Surveys,  X = 7.6  inches, 


Figure  G-5:  Static  Pressure  Surveys,  X = 5.1  inches, 


Figure  G-6:  Static  Pressure  Surveys,  X = 7.6  inches 
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7.599  INCHES 


Static  Pressure  Surveys,  X = 7.6  inches. 
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Figure  G-10:  Static  Pressure  Surveys,  X = 5.1  inches 


Figure  H-2:  Total  Temperature  Surveys,  X = 7.6  inches. 
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Figure  H-3:  Total  Temperature  Surveys,  X = 7.6  inches 
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Total  Temperature  Surveys,  X = 7.6  inches, 
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Figure  H-5:  Total  Temperature  Surveys,  X = 5.1  inches, 
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Figure  H-6:  Total  Temperature  Surveys,  X = 7.6  inches. 
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Figure  H-7:  Total  Temperature  Surveys, 
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Figure  H-8:  Total  Temperature  Surveys,  X = 7.6  inches, 
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7.60  INCHES  , ALFA  G = 10  DEGREES 


ALFA  G = 4 DEGREES 


7.60  INCHES  L ALFA  G = 10  DEGREES 


TOTAL  TEMPERATURE  RATIO  : TT/TSTAG 
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